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S ■ Abstract 

' We perform a QCD sum rule study of the open-charmed Ds(2317) as a four-quark state. Using 

Q ■ the diquark-antidiquark picture for the four-quark state, we consider four possible interpolat- 

. ing fields for Ds(2317), namely, scalar-scalar, pseudoscalar-pseudoscalar, vector- vector, and axial- 

, vector-axial- vector types. We test all four currents by constructing four separate sum rules. The 

\ sum rule with the scalar-scalar current gives a stable value for the Dg mass which qualitatively 

' agrees with the experimental value, and the result is not sensitive to the continuum threshold. The 
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vector-vector sum rule also gives a stable result with small sensitivity to the continuum threshold 



■ and the extracted mass is somewhat lower than the scalar-scalar current value. On the other hand, 



the two sum rules in the pseudoscalar and axial- vector channels are found to yield the mass highly 



\ sensitive to the continuum threshold, which implies that a four-quark state with the combination 

^ \ of pseudoscalar-pseudoscalar or axial-vector-axial-vector type would be disfavored. These results 

would indicate that Z)s(2317) is a bound state of scalar-diquark and scalar-antidiquark and/or 



vector-diquark and vector-antidiquark. 
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I. INTRODUCTION 



Recently new resonances with open- and hidden-charmed states have been reported by 
various experiments. The reported open-charmed states include £)(2308) [1], 1^^.(2317) [2], 
^,(2460) [3, 4], and ^^(2632) [5], while X(3872) [6] and y(4260) [7] have been reported as 
possible hidden-charmed states. This series of new resonances in charm sector opens a new 
challenge for heavy quark system and a systematic analysis is required to understand their 
internal structure and their properties. 

Among the open-charmed states, Ds{2317), which is believed to have the quantum num- 
ber = [8], is particularly interesting. It might be interpreted as a two-quark state of 
cs [9-12], a four-quark state [13-16], a molecular state of DK [17], or a Ds(1968)7r atom [18]. 
The exotic possibility like a four-quark state is quite intriguing as there has been no direct 
and definitive observation of four-quark states so far.^ In this picture, Ds(2317) may be con- 
sidered as a diquark-antidiquark bound state as the diquark picture has been quite useful 
in describing baryon spectroscopy, static properties, and decay mechanisms. (See Ref. [22] 
for a general review on the diquark.) 

Indeed, Bracco et al. [16] recently performed the QCD sum rule calculation of ^5(2317) 
using the current of scalar-diquark-scalar-antidiquark. Their sum rules give the mass which 
agrees with the experimental value. This suggests that the Z}s(2317) may be regarded as a 
four-quark state, especially a bound state of scalar- diquark and scalar-antidiquark. However, 
there can be other possible currents for the diquark-antidiquark system. For example, if the 
light [ud\ diquark is isoscalar, it can have three possible spin quantum numbers: scalar, 
pseudoscalar, and vector [23, 24]. To form a four-quark state with = O"*" in a simple 
approach, therefore, one can combine two diquarks with such quantum numbers, which yields 
three possible choices for the Dg currents, namely, scalar-scalar, pseudoscalar-pseudoscalar, 
and vector- vector. Since the diquarks in the four-quark state of our concern are isodoublet, 
other choices like axial-vector-axial-vector combination may not be excluded. In the light 
quark sector, the axial- vector current gives isovector diquark [24]. 

Of course, in the constituent quark picture the scalar-diquark may be favored over the 
pseudoscalar diquark, because in the pseudoscalar channel the upper component of one Dirac 
spinor is connected only with the lower component of the other spinor. Such combinations 
should vanish in the nonrclativistic limit. However. QCD sum rules [25] deal with current 
quarks and the picture based only on the constituent quarks should be carefully examined. 
In particular, the standard nucleon current [26] involves a linear combination of two types 
of diquarks, scalar and pseudoscalar, and the both contribute to the nucleon sum rule with 
a similar weight. Therefore, a further test with the pseudoscalar-pseudoscalar current is 
required for the Ds sum rule. The other four-quark currents with vector-vector and axial- 
vector-axial-vector types are interesting as they cannot be ruled out from the constituent 
quark picture. Thus, dynamical calculations are needed to test these currents, which would 
be important to understand the possible internal structure of Dg meson when viewed as a 
tetra-quark state. This will be eventually helpful to investigate the Dg properties in lattice 
QCD calculations [27]. 

In this work, we test all the four currents for Ds(2317) in QCD sum rules. For this 



The four-quark state picture of the hght scalar meson nonet of (cr, k, ao, /o) and exotic isotensor mesons 
were discussed, e.g., in Refs. [19, 20]. A possible test for the structure of the light scalar mesons is 
suggested in Ref. [21]. 
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purpose, we first improve the previous sum rule calculation of Ref. [16] by including higher 
order terms in the operator product expansion (OPE) that might be non-negligible for sum- 
rule predictions. We then construct three more sum rules with pseudoscalar-pseudoscalar, 
vector-vector, and axial-vector-axial-vector currents. By scrutinizing those sum rules, we 
hope to eliminate certain diquark combinations as the main composition of the tetraquark 
Dsi which may give a clue for the internal structure of Dg meson. 



II. QCD SUM RULES FOR r>a(2317) 

The quantum numbers of £'5(2317) are believed to be I{J^) = 0(0+) [8]. The four-quark 
current with these quantum numbers can be constructed by combining the diquark [cit] and 
the antidiquark [su\. (For the QCD sum rule calculation for the normal Z}s(1968) meson, 
see, e.g., Ref. [28].) The isoscalar diquark in the light {u^d) sector is restricted to the three 
types: scalar, pseudoscalar, and vector [23, 24]. Thus, its straightforward extension leads to 
the corresponding three types for [cu] and [su\. In addition, the axial- vector diquark may 
be another possibility. By combining the diquark and antidiquark with the same type to 
make the — 0"*" state, we have four possible choices for the four-quark current, 

-^fe = eade^afgC^^kUeSffkU^g + {u ^ d) , (l) 

where /c = 1, . . . , 4 and f ^ = 7°r^7°. The color indices are represented by the subscripts and 
they are chosen so that the diquark in color space belongs to 3c and the antidiquark to 3c. 
The Dirac matrices between the quarks can be Fi = C (pseudoscalar), F2 = C75 (scalar), 
F3 = (^757^ (vector), and F4 — C^^ (axial- vector), where C is the charge conjugation 
operator. For the vector and axial-vector cases, the two vector indices must be contracted 
to form a scalar state. Note that C, C'-f^'-fi, and C70 are off-diagonal matrices and therefore 
they connect the upper and lower components of participating Dirac spinors. The other 
matrices, C75, C7570, and Cji are diagonal so only the upper (or lower) components of the 
Dirac spinors can be connected to each other. Note that, to form an isoscalar four-quark 
current, the {u d) term is added in Eq. (1). 

In constructing sum rules, we consider the following correlation function, 

Uk{p) = i J d^xe'P^iO I T[Jk{x)jl{0)] I 0). (2) 

The time-ordering of quark fields is evaluated by the OPE. This OPE is then matched with 
the hadronic expression of the correlation function via the Borel- weighted sum rule, 

dp^e-P"/^"-Im [nP^-(p2) - nOPE(p2)] = 0, (3) 



J m 



where M denotes the Borel mass and ttIc the charm quark mass. The phenomenological 
side W^^^{p'^) contains the contribution from the low- lying resonance of our concern as well 
as higher resonances or multi-meson continuum states. The contributions other than the 
low-lying resonance have been subtracted according to the QCD duality assumption, which 
introduces the continuum threshold Sq. With this relation, we can extract the Ds mass. 



rrijy , from 



A|2e— = / dp^e-f'/^ -Im [nO™(p2)] , (4) 
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where A represents the coupling strength of the interpolating field to the physical Dg state. 
Specifically, we take a derivative with respect to 1/ and divide the resulting equation by 
Eq. (4) to get the sum rule for m|,^ . Depending on the interpolating fields characterized by 
the index k, we can construct four separate sum rules. To get a sensible result, one has to 
check that the right-hand side of Eq. (4) is positive as constrained by the left-hand side. 

The OPE calculation can be done straightforwardly using the same technique developed 
in Ref . [29] . Namely, we use the momentum-space expression for the charm-quark propagator 
to keep the charm quark mass finite. For the light-quark part, we calculate in the coordinate- 
space, which is then Fourier-transformed to the momentum space in the £)-dimension. The 
resulting light-quark part is combined with the charm-quark part before dimensionally reg- 
ularized at D = 4. The OPE for the scalar-scalar correlator H^^^ is obtained by summing 
up the following terms. 
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(5) 



Here L = u{l — u)p^ — um^ and A = 1 — ml/p^. Compared with the OPE of Ref. [16], 
the terms containing ml and ms{ss) are new, but we found that their contribution is not 
important. In addition, we have higher order OPE containing quark condensate like {qqY 
and {qq) {qGq) , which can constitute nontrivial contributions in stabilizing the Borel curves. 
Using the relation (qD'^q) — {qggaGq) /2, one can combine the two terms at dimension 5. We 
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FIG. 1: Diagrams that contribute to the gluon condensate. The sohd Hnes denote quark propaga- 
tors with specified flavor and the dashed hnes are gluons. When the gluon condensate is formed 
by taking one gluon from one diquark and the other gluon from the other diquark, we have four 
possibilities, (d) shows one of such possibilities. 

have written these terms separately because they are differently related to the corresponding 
terms in the vector and axial-vector channels. Figure 1 shows the diagrams that contribute 
to the gluon condensate. In the equations above, we separate the gluon contribution into 
two terms, one for Figs. 1(a) and (b), and the other for Fig. 1(c). Note that the contribution 
from Fig. 1(d) vanishes in the scalar-scalar channel but it is non-zero in the vector- vector 
channel. We did not calculate the diagrams where two gluons emitted from a light quark 
propagator as they should be parts of the quark condensate. 

For the correlators with the other currents, Hi (pseudoscalar), (vector), and 114 
(axial- vector), most of the OPE can be obtained from 112 by careful inspection. For the 
pseudoscalar correlator, we find 
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while we obtain for vector and axial-vector correlators 
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FIG. 2: (Color online) The Borel curves for the Dg mass in the scalar and pseudoscalar channels. 
The solid curves are for the 112 (scalar) sum rule. The dashed curves are for the Hi (pseudoscalar) 
sum rule. In each case, the thick line is for ^/SQ = 2.7 GeV and the thin line is for \/So = 2.4 
GeV. We observe that for the Hi sum rule the extracted mass is highly sensitive to the continuum 
threshold. But for the 112 sum rule the Borel curves are stable and not so sensitive to the continuum 
threshold. 
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In the above equations, Eg 4 



n 



4n 
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'"^ ^ and H^^f'^^ are calculated directly and their relations to those of 



the scalar OPE look different from the others. In addition, the {qq) {qGq) term in the vector 
and axial-vector correlator is not simply related to that of the scalar correlator. Explicitly, 
it reads 
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III. RESULTS AND DISCUSSION 



Given in Fig. 2 are the Borel curves for the mass m^^ in the scalar-scalar and pseudoscalar- 
pseudoscalar cases. They are obtained by using the following standard values for the QCD 
parameters as used in Ref. [29], 



G^) = (0.33 GeV)^ 



rric = 1.26 GeV, 

{qq) = -(0.23 GeV)^ 

{qgsCrGq) = ml{qq), 



= 0.15 GeV, 

{ss) = 0.8 {qq), 
{sgsCrGs) = 0.8 {qgsCrGq) 



(9) 
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FIG. 3: (Color online) The Borel curves for the Dg mass in the vector and axial-vector channels. 
The solid curves are for the Us (vector) sum rule and the dashed curves are for the 114 (axial- vector) 
sum rule. The thick lines are for ^/Sq = 2.7 GeV and the thin lines for ^/Sq = 2.4 GeV. 



The parameter ttiq denotes the quark virtuahty which is normally taken to be 0.8 GeV^ [30]. 
To show the sensitivity to the continuum threshold, we vary ^/So from 2.4 GeV (thin lines) 
to a somewhat larger value, 2.7 GeV (thick lines). The lower bound, 2.4 GeV, is slightly 
above the KD threshold. The 112 sum rule shown by the two solid lines clearly exhibits 
the Borel stability yielding m^^ around 2.4-2.5 GeV depending on the continuum. The two 
solid lines (thin and thick lines) are quite close to each other indicating that this sum rule is 
insensitive to the continuum threshold. Our result is not so sensitive to the uncertainties of 
the most QCD parameters given above except for the parameter associated with the quark- 
gluon mixed condensate, rriQ. If somewhat lower value = 0.6 GeV^ is used in this sum 
rule, the extracted mass with ^/So = 2.4 GeV is found to be around 2.27 GeV. In fact, the 
quark-gluon mixed condensate plays a crucial role in making the right-hand side of Eq. (4) 
be positive as constrained by the left-hand side. We found that, without the quark-gluon 
condensate, this constraint is not satisfied in this particular sum rule. Our stable result may 
imply that the scalar-scalar current couples strongly to the low-lying pole while its couplings 
to higher resonances or continuum states are suppressed. Therefore, the scalar-scalar current 
might be an optimal candidate for the = O"*" Dg resonance. 

On the other hand, the Borel curves for the Hi sum rule given by the dashed lines in 
Fig. 2 show strong sensitivity to the continuum threshold. The thick dashed line {\/So = 2.7 
GeV) is substantially different from the thin dashed line (a/Sq = 2.4 GeV). Also, the strong 
variation of the curves with respect to the Borel mass clearly shows the Borel instability. 
From this sum rule, therefore, it is hard to extract any stable value for the Dg mass. In this 
case, the important contribution coming from the quark-gluon mixed condensate changes 
its sign from that of 112 as can be seen from Eq. (6), which makes the Borel curves unstable. 
Thus, the quark-gluon mixed condensate contributes differently in the pseudoscalar sum 
rule. Note that, when y/So = 2.4 GeV, there is no extracted mass for M > 2.82 GeV as it 
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FIG. 4: The Dg mass as a function of the charm quark mass rric with the scalar-scalar (112) current 
and with the vector- vector (lis) current at the Borel mass M = 2.4 GeV. The shaded area for 112 
current is obtained by varying rriQ from 0.7 GeV^ to 0.9 GeV^, while that for II^ current by varying 
(qq) = —(0.23 GeV)'^ ~ —(0.25 GeV)^. Other parameters arc fixed as given in Eq. (9) in each case 
with = 2.5 GeV. The dashed line denotes the experimental mass of Ds(2317). 

becomes imaginary. One way to understand the success (failure) of the scalar (pseudoscalar) 
sum rule may be the dominance of the nonrelativistic configuration of the current for the 
Dg sum rules. Namely, the scalar current survives in the nonrelativistic limit while the 

pseudoscalar current does not. 

Figure 3 shows the results from the Ho, and 114 sum rules. Wc observe that for the II3 
sum rule the extracted mass is slightly sensitive to the continuum threshold. As we change 



y/So from 2.4 GeV to 2.7 GeV, the extracted mass is in the range of 1.92 to 2.2 GeV, which 

is somewhat smaller than the experimental Ds mass and the value obtained with the 112 
sum rule. But the extracted mass can be larger when a slightly larger continuum threshold 
is used. Furthermore, the Borel curves are quite stable with respect to the Borel mass and 
thus it can be another possible choice for the Dg current. In the 114 sum rule, however, the 
Borel curves are quite sensitive to the continuum threshold as one can see from the dashed 
lines in Fig. 3. The result from the 114 sum rule with ^/Sq = 2.4 GeV is in fact unphysical 
since the OPE is found to violate the positivity constraint provided by Eq. (4). Thus, the 
current with axial-vector-axial-vector composition is ruled out from the possible choices 
for the Dg current. If this result should be related to the dominance of the nonrelativistic 
configuration as conjectured in the scalar and pseudoscalar channels, then our results in 
the vector and axial-vector channels can be understood from the dominance of the time 
component of diquarks over their spatial components. 

As we have discussed above, the obtained masses are mostly dependent of the continuum 
threshold 5*0. In order to estimate the uncertainties arising from the other parameters we 
plot the obtained Dg mass as a function of the charm quark mass rric in Fig. 4. Among the 
parameters given in Eq. (9), we found that the Dg mass with 112 current is most sensitive 
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to the quark virtuality itLq, while the mass obtained with lis is to the quark condensate 
(qq). The shaded area in Fig. 4 are, therefore, obtained with ttIq = 0.7 GeV^ ~ 0.9 GeV^ 
for 112 and with (qq) = —(0.23 GeV)^ ~ —(0.25 GeV)^ for U^. As the results are much less 
sensitive to the other parameters, they are fixed as given in Eq. (9) and \/^ — 2.5 GeV 
is used. Within the range of the parameters considered in this calculation. Fig. 4 shows 
that the uncertainties are not large and actually they are at the order of 5% for a given 
continuum threshold y/So- 




In summary, we have performed QCD sum rule calculations for ^5(2317) using four 
different tetra-quark currents with the final quantum numbers — 0'^. The sum rule using 
the scalar-scalar current gives the stable Borel curves with least sensitivity to the continuum 
threshold. This slightly improves the previous calculation by Bracco et al. [16] and supports 
the picture of 1^3(23 17) as a bound state of scalar-diquark and scalar-antidiquark. Also the 
vector- vector sum rule was found to give a stable result for the Dg mass with slightly lower 
value. This implies that the vector-vector combination could be another possible choice and 
cannot be simply ruled out for the Dg current. Of course, the possibility that the physical 
Ds state may be a mixture of the two configurations cannot be excluded. The other two 
sum rules using the pseudoscalar-pseudoscalar and axial-vector-axial-vector currents yield 
unstable results with strong sensitivity to the continuum threshold. Therefore, they are 
unfavored as the main component of the Ds(2317) meson. Although we could not rule out 
the other models, such as quark-antiquark description,^ for the Ds meson by this study, our 
results indicate the possible four-quark structure with the scalar-scalar and/or vector- vector 
diquarks of the Dg meson. 
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